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THE  FLUID  MECHANICS  OF  VORTEX 
CUTTING  BY  A  BLADE 


1.  Statement  of  the  Problem 


The  problem  under  consideration  in  this  study  involves  the  impingement  of  a  vortex 
onto  a  blade,  where  the  blade  travels  relative  to  the  vortex  in  a  direction  normal  to  the 
vortex  axis,  as  shown  in  figure  1.  This  problem  will  hereinafter  be  referred  to  as  normal 
blade- vortex  interaction  (BVI).  The  primary  application  of  this  research  is  in  the  area  of 
helicopter  interactional  aerodynamics.  During  the  course  of  the  research,  the  scope  of  the 
study  has  evolved  such  that  we  interpret  the  term  "blade",  in  a  loose  manner,  to  be  any 
body  which  is  very  long  in  one  direction,  oriented  approximately  normal  to  the  vortex 
axis  (see  figure  1).  The  research  would  thus  apply  to  problems  involving  interaction  of 
helicopter  main  rotor  vortices  with  the  tail  or  empennage  sections  of  the  helicopter,  which 
is  a  common  occurrence  in  low  speed  flight,  or  (for  some  designs)  with  the  tail  rotor 
blades  (Leverton  et  al.,  1977).  The  research  would  also  apply  to  problems  such  as 
ingestion  of  atmospheric  turbulence  or  the  helicopter  "ground  vortex"  into  the  main  rotor. 

The  results  of  this  project  have  been  reported  in  two  parts,  which  respectively 
consider  the  nature  of  the  vortex  response  to  the  blade  as  a  function  of  the  two 
dimensionless  parameters,  T/a  and  Y HiiaU .  The  first  of  these  parameters  represents 
the  ratio  of  the  blade  "thickness",  T,  to  the  ambient  vortex  core  radius,  a,  while  the 
second  parameter  represents  the  ratio  of  the  maximum  nominal  vortex  swirl  velocity, 
r/27ta,  to  the  blade  speed,  U,  relative  to  the  vortex  axis.  The  effects  of  these 
parameters  on  the  vortex  flow  is  examined  using  a  combination  of  experiment,  numerical 
computation  and  analytical  modeUng. 

The  first  part  of  the  study,  which  has  previously  been  reported  to  ARO  (Marshall, 
1993),  concentrates  primarily  on  the  effect  of  T  /  a  on  the  vortex  response.  It  is  found 
that  this  parameter  is  primarily  responsible  for  determining  the  amount  that  a  vortex  will 
bend,  due  to  the  inviscid  interaction,  as  it  comes  close  to  the  blade.  Two  basic  regimes  of 
interaction  were  identified,  one  for  small  T/a,  in  which  the  vortex  remains 
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approximately  unbent  as  comes  close  to  and  is  eventually  cut  by  the  blade,  and  one  for 
large  Tic,  in  which  the  vortex  undergoes  a  large  displacement  and  wraps  about  the 
blade  leading  edge. 


The  second  part  of  the  research,  reported  here,  is  concerned  with  the  effect  of  the 
second  parameter,  F  /  liicU ,  on  the  vortex  response.  In  the  context  of  an  inviscid  flow 
with  small  Tic,  so  that  the  blade  penetrates  into  the  vortex,  the  combination 
{T !  c){Y  UncU)  is  found  to  play  a  central  role  in  determining  whether  the  vortex 
deformation  is  due  primarily  to  the  vortex-induced  flow  or  to  the  ambient  potential  flow 
past  the  blade.  In  the  latter  case,  approximations  such  as  rapid  distortion  theory  (RDT) 
can  be  used  to  simplify  calculation  of  the  unsteady  force  on  the  blade,  as  is  commonly 
done  in  models  of  sound  generation  from  blade-vortex  interaction  (Howe,  1989,  Amiet  et 
al.,  1990).  For  viscous  flow,  the  parameter  F  /  IwcU  influences  the  time  at  which 
secondary  vorticity  is  shed  from  the  blade  by  the  vortex-induced  flow  and  the  extent  to 
which  the  shed  vorticity  influences  the  primary  vortex  prior  to  impinging  on  the  blade. 
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Fig.  1.  Schematic  of  normal  blade- vortex  interaction,  showing  core  axial  flow  w,  vortex 

core  radius  c,  vortex  circulation  F,  blade  speed  U  and  angle  of  attack  a. 
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2.  Summary  of  the  Most  Important  Results 


The  simplest  model  of  vortex-blade  interaction  is  that  of  an  infinitesimally  thin 
vortex  filament  impinging  on  a  blade  of  thickness  T,  where  the  vortex  filament  wraps 
around  the  blade  and  comes  progressively  closer  to  the  blade  surface  as  time  advances. 
Let  b  denotes  the  minimum  separation  distance  between  the  vortex  filament  and  the  blade 
surface.  As  the  vortex  filament  comes  very  close  to  the  surface  the  blade  surface  pressure 
would  exhibit  a  1  /  singularity  and  the  force  exerted  on  the  blade  by  the  vortex  would 
exhibit  a  1/Z?  singularity,  thus  becoming  infinite  as  Z?  0 .  This  singularity  is  clearly  non¬ 
physical;  it  is  an  artifact  of  the  assumption  in  the  model  that  all  of  the  vorticity  is 
concentrated  on  a  curve. 

There  are  several  ways  in  which  the  force  singularity  might  be  resolved.  For 
instance,  one  might  assume  that  the  vortex  core  remains  circular  but  has  a  finite  radius 
which  is  uniform  along  the  vortex  (Affes  and  Conlisk,  1993).  One  might  alternatively 
make  allowance  for  the  fact  that  the  part  of  the  vortex  near  the  blade  is  stretched  more 
rapidly  than  the  part  far  from  the  blade  by  allowing  the  core  radius  of  the  filament  to  vary 
with  position  along  the  vortex  axis,  but  then  it  is  also  necessary  to  account  for  the  strong 
axial  flow  within  the  vortex  that  forms  as  a  result  of  the  gradient  in  core  radius  (Marshall 
and  Y2damanchili,  1994).  One  might  let  the  core  shape  deform  from  a  circular  shape, 
which  would  occur  due  to  interaction  with  the  vortex  image  (e.g.,  Saffman  and  Tanveer, 
1982).  For  cases  where  the  blade  thickness  is  of  the  same  order  as  or  much  thinner  than 
the  vortex  core  radius,  the  blade  is  observed  to  penetrate  into  the  vortex  core,  in  which 
case  it  is  necessary  to  account  for  the  very  strong  vorticity  deformation  occurring  close  to 
the  leading  edge  of  the  blade  (Marshall  and  Grant,  1995).  Viscous  effects  can  also 
substantially  modify  the  vortex-blade  interaction  if  the  vorticity  ejected  from  the  blade 
boundary  layer  has  sufficient  time  to  interact  strongly  with  the  primary  vortex  (Kim  and 
Komerath,  1995). 

The  important  point  to  be  made  in  the  above  discussion  is  that  the  maximum  force 
induced  by  the  vortex  on  the  blade  is  strongly  effected  by  details  of  the  vortex  structure 
and  deformation  during  impact  of  the  vortex  on  the  blade.  Research  in  this  area, 
particularly  with  respect  to  acoustic  aspects  of  vortex-blade  interaction,  typically  employs 
some  model  of  the  vortex  flow  (such  as  that  of  a  "frozen"  vortex,  inviscid  flow,  rapid 
distortion  theory,  circular  core,  etc.),  but  it  is  not  at  all  clear  that  the  choice  of  model  does 
not  have  a  significant  effect  on  the  predicted  maximum  blade  force  and  sound  generation. 
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On  the  other  hand,  direct  computation  of  vortex  impact  on  a  blade  in  three  dimensions 
from  the  Navier-Stokes  equations  is  quite  difficult  and  is  not  likely  to  be  compatible  with 
time  requirements  of  the  helicopter  design  process.  Some  form  of  modeling  is  clearly 
necessary  if  research  in  this  area  is  to  have  any  impact  on  helicopter  design. 

The  approach  taken  in  the  present  study  is  to  first  identify  the  physical 
characteristics  of  the  flow  during  vortex  impact  on  the  blade  in  different  parameter 
regimes,  and  then  to  evaluate  possible  computational  and  analytical  models  which  can  be 
used  to  describe  aspects  of  the  flow  field.  Section  2a  gives  an  overview  of  the  physical 
characteristics  of  the  flow  field  during  normal  BVI  as  a  function  of  the  two  parameters 
T !  o  and  F  UncsU .  A  complementary  summary  was  given  by  Kim  and  Komerath 
(1995),  which  describes  various  phenomena  associated  with  rotor  wake-cylinder 
interaction  but  does  not  span  parameter  space  of  F  /  a  and  F  /  InaU .  Also,  Sheridan 
and  Smith  (1979)  discuss  the  need  for  research  in  interactional  aerodynamics  in  terms  of 
helicopter  flight  performance.  Section  2b  considers  the  case  in  which  the  blade 
penetrates  into  the  vortex  core,  and  provides  an  assessment  of  the  RDT  approximation  for 
this  problem,  which  forms  the  basis  of  many  acoustic  analyses  for  normal  vortex-blade 
interaction  (Howe,  1989,  Amiet  et  al.,  1990).  Conclusions  and  suggestions  for  future 
work  are  given  in  section  2c. 

2a.)  Regimes  of  normal  blade-vortex  interaction 

Experiments  have  been  performed  in  which  a  thin  blade  or  a  circular  cylinder  is 
moved  through  a  nominally  straight  vortex  with  a  non-zero  axial  flow.  The  vortex  was 
formed  using  the  vortex  generator  device  described  by  Krishnamoorthy  and  Marshall 
(1994),  using  water  as  the  working  fluid.  The  diameter  of  the  circular  cylinder,  the 
towing  velocity  of  the  blade  or  cylinder  and  the  flow  rate  through  the  vortex  generator 
device  could  all  be  varied,  such  that  the  parameters  T  /  a  and  F  /  2kgU  could  be 
adjusted  from  0.7  -  17.6  and  from  3.3  -  195,  respectively.  The  axial  velocity  w  in  the 
vortex  for  all  cases  considered  was  such  that  the  ratio  2 new  /  F  was  equal  to  about  0.15. 
The  angle  of  attack,  a,  of  the  blade  could  also  be  adjusted  from  -15°  to  -1-15°. 
Experimental  results  were  obtained  primarily  by  marking  the  intake  vortex  with  red  dye 
and  the  blade  boundary  layer  by  blue  dye,  and  then  photographing  the  resulting  flow 
evolution.  Measurements  of  vortex  circulation  and  vorticity  distribution  were  obtained 
both  by  photographing  the  motion  of  small  (1-2  mm)  immiscible  dye  globules  immersed 
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in  the  boundary  layer  (Krishnamoorthy  and  Marshall,  1994)  and  by  particle-image 
velocimetry. 

(i)  Regimes  of  vortex-body  interaction 

Three  qualitatively  different  regimes  of  vortex-blade  (or  cylinder)  interaction  were 
identified.  The  first  regime  occurs  for  cases  in  which  the  parameter  T  l2noU  is  large 
(greater  than  about  20-30),  so  that  the  swirl  velocity  of  the  vortex  is  much  larger  than  the 
convection  speed  of  the  vortex  axis  relative  to  the  blade.  In  this  strong-vortex  regime, 
vorticity  is  ejected  from  the  blade  boundary  layer  while  the  vortex  is  still  quite  distant 
from  the  blade  (5-10  core  diameters  or  more).  A  sequence  of  pictures  showing  the 
evolution  of  the  ejected  vorticity  for  a  case  with  Tic  =  17.6  and  T  UtzcU  =  40  is 
given  in  figures  2a-d.  The  ejected  vorticity  rolls  up  into  a  vortex  filament  which  is 
attached  to  the  blade  boundary  layer  at  two  points.  This  vortex  filament  wraps  about  the 
primary  vortex,  and  the  induced  velocity  from  the  secondary  vortex  loops  cause  formation 
of  both  axial  waves  (associated  with  variation  of  the  core  area)  and  helical  waves 
(associated  with  displacement  of  the  core  axis)  on  the  primary  vortex.  The  wavelength  of 
these  waves  on  the  primary  vortex  is  of  the  same  order  as  the  distance  between  secondary 
vortex  loops.  These  waves  continue  to  grow  in  amplitude  and  decrease  in  wavelength, 
eventually  resulting  in  complete  disruption  of  the  primary  vortex,  which  may  occur  (for 
large  T  UncU)  while  the  primary  vortex  is  still  separated  from  the  blade  by  several  core 
diameters. 

The  second  regime  occurs  for  weak  vortices  (T  /  lucU  less  than  about  20)  with 
thin  blades  {Tic  less  than  about  2).  In  this  regime,  the  vorticity  ejected  from  the  blade 
boundary  layer  has  little  time  to  influence  the  primary  vortex,  and  the  blade  penetrates 
into  the  vortex  core  with  no  significant  deflection  of  the  primary  vortex  axis.  Once  the 
blade  penetration  into  the  core  occurs,  computations  for  inviscid  flow  (see  section  2b) 
indicate  that  the  vorticity  will  accumulate  in  a  thin  "ribbon"  of  high-vorticity  fluid  which 
wraps  about  the  blade  leading  edge.  This  vortex  ribbon  eventually  decays  due  to  viscous 
interaction  with  the  underlying  boundary  layer  on  the  blade,  leading  to  connection  of 
vorticity  lines  from  the  primary  vortex  with  those  in  the  blade  boundary  layer. 
Experiments  with  two  color  dyes  indicate  a  strong  entrainment  of  boundary  layer  fluid 
into  the  primary  vortex  core  both  above  and  below  the  blade. 
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Pictures  showing  secondary  vortex  wrapping  about  the  primary  vortex  in  the 
strong-vortex  regime  (regime  1).  (a)  Vorticity  is  ejected  from  blade  boundary 
layer  and  pulled  away  from  the  blade  by  the  primary  vortex,  (b)  Secondary 
vortex  filament  completes  one  wrap  about  the  primary  vortex. 
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Fig.  2.  Pictures  showing  secondary  vortex  wrapping  about  the  primary  vortex  in  the 
strong- vortex  regime  (regime  1).  (c)  Secondary  vortex  filament  completes  one 
and  one-half  wrap  about  primary  vortex,  and  waves  begin  to  form  on  the 
primary  vortex,  (d)  Many  waves  form  on  primary  vortex  due  to  wrapped 
secondary  vorticity,  and  the  vortex  disrupts  completely  shortly  thereafter. 
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The  third  regime  occurs  for  weak  vortices  (T  /  IncsU  less  than  about  10)  but  thick 
bodies  (T/a  greater  than  about  5).  In  this  regime,  the  vorticity  ejected  from  the 
boundary  layer  again  has  little  time  to  influence  the  primary  vortex  before  it  impacts  on 
the  body.  However,  for  thick  bodies  (large  T  /  c),  the  vortex  is  observed  both  in 
experiments  and  in  numerical  computations  (Affes  and  Conlisk,  1993,  Marshall  and 
Yalamanchili,  1994)  to  bend  about  the  blade  surface,  as  shown  in  figure  3  with 
T/g  =  \0.  For  very  small  values  of  T  /InoU,  substantial  thinning  of  the  vortex  core 
occurs  due  to  stretching  by  the  blade  (Marshall  and  Yalamanchili,  1994).  After  the 
vortex  bends  about  the  blade  a  certain  amount  (in  a  manner  qualitatively  similar  to  that 
predicted  by  vortex  filament  calculations),  small  helical  waves  develop  on  the  vortex 
axis.  These  helical  waves  grow  and  eventually  impinge  on  the  boundary  layer  of  the 
body,  at  which  point  the  vortex  rapidly  disrupts  into  small  scale  turbulence.  These  helical 
waves,  which  lead  to  breakup  of  the  vortex,  may  be  caused  either  by  inviscid  instability  of 
the  vortex  with  its  image  over  the  body  surface  (Pedrizzetti,  1992)  or  by  interaction  with 
ejected  vorticity  from  the  body  boundary  layer. 

In  all  three  regimes  of  vortex-body  interaction,  a  vortex  breakdown  forms  on  the 
vortex  during  impact  with  the  blade  and  propagates  up  the  vortex  axis,  opposite  to  the 
direction  of  the  ambient  axial  flow  in  the  core.  The  formation  of  the  propagating  vortex 
breakdown  following  cutting  of  the  vortex  is  consistent  with  predictions  of  formation  of  a 
propagating  vortex  "shock"  (i.e.,  a  discontinuity  in  vortex  core  area)  obtained  using  an 
extended  vortex  filament  model  which  allows  for  variation  in  vortex  core  area  (Marshall, 
1994),  as  illustrated  in  figure  4.  The  measured  propagation  speed  of  the  traveling  vortex 
breakdown  also  agrees  well  with  that  predicted  for  the  vortex  "shock"  in  the  filament 
model  (Krishnamoorthy  and  Marshall,  1994).  The  traveling  vortex  breakdown  was 
observed  to  be  either  of  a  spiral-  or  bubble-form  (shown  in  figure  5a  and  b),  whieh  occurs 
for  low  and  high  values  of  T/2naU,  respectively.  The 

dependence  of  breakdown  form  on  T  /  InoU ,  observed  in  the  experiments,  might  be  due 
to  the  fact  that  the  breakdown  develops  when  the  vortex  is  still  detached  from  the  blade 
for  large  values  of  T  /  27taC/,  while  for  small  values  of  F  /  27cal7  it  develops  when  the 
blade  penetrates  into  the  vortex  core. 

(»■)  Vortex-induced  forces  on  the  blade 

The  vortex-induced  forces  on  the  blade  during  normal  BVI  result  from  three 
qualitatively  different  effects:  (1)  the  reduction  in  surface  pressure  due  to  swirling  motion 
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Fig.  5.  Pictures  of  vortex  breakdowns  of  the  bubble  and  spiral  types,  both  generated  by 
cutting  the  vortex  with  a  thin  blade  at  zero  angle  of  attack  and  propagating 
upward  on  the  vortex  axis.  A  bubble-type  breakdown,  which  spirals  into  a 
double-helix  downstream,  is  shown  in  (a)  for  the  case  27tGw  /  F  =  —0.16  and 
r  /  27ia<7  =  28.  A  spiral-type  breakdown  is  shown  in  (b)  for  the  case 
2k(Sw  /  r  =  —0.26  and  F  /  litoU  =  5.3. 


12 


of  the  vortex,  (2)  blockage  of  the  vortex  axial  flow  as  the  blade  cuts  the  vortex  core,  and 
(3)  vorticity  accumulation  on  the  blade  surface  when  the  blade  penetrates  into  the  vortex 
core.  The  first  and  the  third  of  these  forces  are  oriented  in  the  direction  of  relative  motion 
between  the  blade  and  the  vortex  axis,  and  the  second  force  is  oriented  in  a  direction 
tangent  to  the  vortex  axis. 

The  first  of  these  forces  can  be  accurately  estimated  using  various  vortex  filament 
theories  (Affes  and  Conlisk,  1993,  Affes  et  al.,  1993,  Marshall  and  Yalamanchili,  1994), 
provided  that  the  vortex  is  separated  from  the  blade  surface  by  several  core  diameters. 
When  the  separation  distance  between  the  vortex  and  the  blade  is  on  the  order  of  a  core 
diameter  or  less,  deformation  of  the  core  and  ejection  of  vorticity  from  the  surface  may 
significantly  influence  the  blade  force. 

A  model  for  the  second  force,  due  to  blockage  of  axial  flow,  is  given  by  Marshall 
(1994)  and  numerical  calculations  for  this  problem  assuming  axisymmetric  inviscid  flow 
have  been  performed  by  Lee  et  al.  (1995),  but  the  forces  predicted  by  these  studies  have 
not  yet  been  verified  experimentally.  While  the  predictions  of  the  model  of  Marshall 
(1994)  concerning  formation  and  propagation  of  vortex  breakdown  on  one  side  of  the 
vortex  after  it  has  been  cut  by  the  blade  agree  well  with  experiments,  this  breakdown  is 
observed  in  experiments  to  have  a  form  which  is  sometimes  far  from  axisymmetric  (as 
shown  in  figure  5b).  The  assumption  of  axisymmetry  for  prediction  of  the  axial  flow 
blockage  force  thus  requires  further  examination. 

The  third  force  mentioned  above,  which  we  will  refer  to  simply  as  the  blade 
penetration  force,  occurs  only  for  the  second  regime  of  blade-vortex  interaction,  for 
which  T  /  a<  0(1)  and  T  /  2noU  is  less  than  about  10.  Since  T  /  litaU  is  relatively 
small  in  this  regime,  it  is  plausible  that  secondary  vorticity  ejection  does  not  play  a  major 
role,  so  that  the  problem  can  be  considered  (at  least  initially)  in  terms  of  inviscid  flow 
theory.  The  solution  for  the  penetration  force  is  greatly  simplified  by  making  use  of  rapid 
distortion  theory  (RDT),  in  which  it  is  assumed  that  the  vorticity  field  is  convected  only 
by  the  ambient  potential  flow  about  the  blade.  The  use  of  RDT  for  calculation  of  blade 
penetration  force  was  examined  by  comparison  to  direction  computation  by  Marshall  and 
Grant  (1995),  and  the  results  of  this  investigation  are  reported  in  the  next  section. 
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2b.)  Computations  of  blade  penetration  into  a  vortex  core 


Numerical  calculations  were  performed  for  penetration  of  a  NACA  series  0016 
airfoil,  with  equal  span  and  chord  lengths,  into  the  core  of  a  vortex  ring.  A  uniform  flow 
with  speed  U  carries  the  vortex  ring  into  the  blade  leading  edge.  All  length  and  time 
scales  are  nondimensionalized  by  the  chord  length  c  and  the  ratio  dU,  respectively.  The 
core  radius,  a,  is  0.1  and  the  ring  radius,  R,  is  0.5.  The  initial  configuration  for  the 
calculations  is  shown  in  figure  6.  A  coordinate  system  is  introduced  such  that  the  z- 
direction  is  along  the  blade  leading  edge,  the  jc-direction  is  along  the  blade  chord  and  the 
y-direction  is  normal  to  the  symmetry  plane  of  the  blade.  The  blade  is  stationary  and  its 
leading  edge  is  coincident  with  the  line  x  =  y  =  0 .  The  center  of  the  ring  is  initially 
located  at  a  distance  of  x  =  -0. 75  upstream  of  the  blade  leading  edge.  Since  the  ring  also 
propagates  in  the  z-direction  under  its  self-induced  velocity  field,  an  initial  spanwise 
displacement  of  the  ring  is  introduced  such  that  the  blade-vortex  collision  will  occur  close 
to  the  blade  center  (z  =  0). 

(0  Numerical  method 

The  vortex  evolution  was  computed  by  direct  solution  of  the  Euler  equation  for 
inviscid  flows  in  vorticity-velocity  space.  The  vorticity  in  the  flow  was  marked  using 
some  number  N  control  points.  All  of  the  calculations  reported  here  were  performed  with 
both  moderate  resolution  {N  =  3500)  and  high  resolution  {N  =  10,500),  and  the  control 
points  were  three  times  as  closely  spaced  on  the  side  of  the  vortex  which  impacted  on  the 
blade.  The  control  points  were  convected  by  the  fluid  flow,  and  the  vorticity  values  on 
the  control  were  evolved  using  the  inviscid  vorticity  transport  equation.  The  velocity 
field  was  obtained  by  solution  of  the  Biot-Savart  integral,  using  Gaussian  blob  functions 
to  interpolate  the  vorticity  field.  Unlike  traditional  vortex  blob  methods,  the  amplitudes 
of  the  Gaussian  blobs  were  fit  to  the  vorticity  values  at  the  control  points  at  every  time 
step  using  a  new  iterative  procedure,  which  is  subjected  to  numerous  validation  tests  by 
Marshall  and  Grant  (1995).  The  no-penetration  boundary  condition  on  the  blade  surface 
was  enforced  using  vorticity  sheet  panels  on  the  surface  (Hess  and  Smith,  1967),  where 
the  current  computations  utilize  900  quadrilateral  panels  to  discretize  the  blade  surface. 
The  blade  surface  pressure  was  solved  using  an  integral  equation  formulation  of  the 
pressure-Poisson  equation  (Uhlman,  1992). 
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Time  advancement  of  the  system  was  performed  using  the  second-order  predictor- 
corrector  method  with  a  time  step  of  0.005.  The  calculations  were  run  for  70  time  steps 
and  were  stopped  at  time  r  =  0.35,  at  which  point  the  blade  has  completely  penetrated 
through  the  ambient  position  of  the  vortex  core.  For  the  moderate  resolution  cases,  runs 
were  performed  in  which  the  time  step  was  allowed  to  vary  during  the  calculations 
between  about  0.002  and  0.0001,  with  no  noticeable  difference  in  the  computed  results. 

The  radius  of  the  vorticity  elements  was  fit  adaptively,  by  searching  for  nearest 
neighbors  in  three  directions  about  each  element,  in  order  to  ensure  sufficient  overlap 
between  elements.  The  nominal  initial  element  radius  size  was  0.04  for  the  moderate 
resolution  cases  and  0.03  for  the  high  resolution  cases.  The  element  radius  decreased  for 
control  points  that  approached  the  blade  surface  to  about  0.03  -  0.025  for  the  moderate 
resolution  cases  and  about  0.015  -  0.008  for  the  high  resolution  cases. 

(ii)  Vorticity  evolution 

In  the  present  computations,  the  value  of  T  /  a  is  fixed  at  1.6,  and  computations 
are  performed  for  three  values  of  the  parameter  F  /  2%aU ,  given  by  0.32,  0.79  and  1.6. 
All  of  these  values  lie  in  the  second  regime  of  vortex-blade  interaction  described  in 
section  2a,  so  we  would  expect  the  blade  to  penetrate  into  the  vortex  core  with  little 
deflection  of  the  vortex  axis.  One  objective  of  these  calculations  is  to  compare  results 
using  the  RDT  approximation  to  those  using  the  full  inviscid  theory.  In  RDT,  the 
vorticity  is  assumed  to  be  deformed  by  only  the  ambient  potential  flow  about  the  blade, 
but  the  full  velocity  field  (including  that  induced  by  the  vortex  and  the  vortex  image) 
must  be  including  in  computing  the  blade  surface  pressure.  In  order  for  RDT  to  be  valid, 
it  is  necessary  for  the  time  scale  t^  ~T !  U  associated  with  deformation  of  the  vorticity 
by  the  ambient  potential  flow  to  be  small  compared  to  the  convection  time  4  =  2;ca^  /  F 

of  the  vortex-induced  motion.  This  condition  implies  that  the  parameter 
T  r 

e  =  tj  I  tc  ~  i^)^2'KGU^  ^ considered  in  the  paper,  with  T /  a  =  1.6 

and  r/27tt/a  =  0.32,  0.79  and  1.6,  correspond  to  8  values  of  0.5,  1.3  and  2.5, 
respectively. 

In  order  to  evaluate  the  ability  of  the  numerical  method  to  resolve  the  vorticity 
deformation  throughout  the  blade  penetration  process,  the  vortex  circulation  is  calculated 
over  a  cross-section  of  the  vortex  core  in  the  plane  y  =  0,  which  coincides  with  the  plane 
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occupied  by  the  blade  chord.  The  circulation  calculation  is  performed  using  two  different 
methods.  For  the  first  method,  the  vortex-induced  velocity  tangent  to  a  circuit 
surrounding  the  core,  with  a  radius  of  three  times  the  core  radius,  is  integrated  over  the 
circuit  to  yield  circulation.  In  the  second  method,  the  vorticity  is  interpolated  using  the 
vortex  blob  representation  onto  a  uniform  grid  on  the  plane  y=0,  with  151  grid  points  on 
each  side  spanning  the  interval  (-0.5,0.05)  in  x  and  (-0.4, 0.4)  in  z,  and  the  component 

of  vorticity  normal  to  this  plane  is  integrated  over  area  to  obtain  the  circulation.  A  plot 
showing  the  results  for  these  two  methods  of  circulation  measurement  is  given  in  figure  7 
for  high  resolution  runs  with  the  three  values  of  T  considered  in  the  paper.  In  this  figure, 
square  symbols  denote  results  using  the  first  method  for  circulation  measurement  and 
circular  symbols  denote  results  using  the  second  method.  The  dashed  lines  correspond  to 
the  nominal  values  of  circulation  prescribed  in  the  initial  plane  for  the  three  cases.  For  all 
times  shown,  except  for  the  last  time  considered  for  the  case  F  =  1  (or  F  /  2nUo  =  1.6), 
the  circulation  obtained  by  either  of  these  methods  differs  from  the  initial  circulation  by 
no  more  than  10%.  In  the  last  point  shown  in  figure  7  for  the  case  F  =  1,  the  circulation 
obtained  by  the  second  method  is  37%  larger  than  the  initial  value,  whereas  that  obtained 
by  the  first  method  is  about  10%  larger  than  the  initial  value.  This  point  occurs  at  a  time 
when  the  penetration  of  the  blade  through  the  ambient  position  of  the  vortex  core  is  87% 
complete.  The  differences  between  the  two  measures  of  circulation  at  this  point  could  be 
indicative  either  of  inadequacy  of  the  element  resolution  during  the  later  stages  of  the 
penetration  (as  will  be  discussed  further  presently)  or  of  significant  non-zero  values  of  the 
vorticity  divergence. 

It  is  found  for  all  three  values  of  F  /  InUa  considered  that  the  vortex  core 
exhibits  relatively  little  response  as  it  approaches  the  blade,  aside  from  a  slight  flattening 
of  the  core,  until  the  point  where  the  blade  leading  edge  actually  begins  to  penetrate  into 
the  vortex  core.  If  the  vortex  were  to  be  convected  with  the  uniform  flow  speed  alone, 
the  vortex  would  begin  penetrating  the  blade  at  a  dimensionless  time  of  t  =  0.15  and 
penetrate  entirely  through  the  blade  by  t  =  0. 35 .  A  plot  showing  the  vorticity  field  for  the 
case  F /  2tzgU  =  1.6  for  a  moderate  resolution  run  at  dimensionless  time  t  =  0.30,  when 
the  penetration  process  is  about  three-quarters  complete,  is  given  in  figure  8.  As  would 
be  expected  from  Helmholtz's  observation  that  vortex  lines  and  material  lines  coincide  for 
an  inviscid  fluid,  the  computed  vortex  lines  deform  about  the  blade  leading  edge  but  are 
not  cut  by  the  blade.  The  highest  vorticity  magnitudes  are  observed  in  the  computations 
to  occur  very  close  to  the  blade  surface  and  to  form  a  'ribbon'  of  vorticity  about  the 
surface  as  the  vortex  is  convected  past  the  blade. 
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circulation 


Fig.  7.  Plot  of  the  variation  with  time  of  the  computed  values  of  circulation  in  the 
y  =  0  plane  as  the  vortex  is  penetrated  by  the  blade  for  the  three  nominal 
circulation  values  considered  (shown  by  a  dashed  line).  Points  denoted  by 
squares  are  computed  using  integration  of  velocity  over  a  contour  about  the 
vortex  core  and  points  denoted  by  circles  are  computed  by  integration  of 
vorticity  over  the  cross-sectional  area  of  the  core. 
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Fig.  8.  Side  view  of  the  vorticity  field  for  the  case  with  F  /  2kgU  =  1.6  at  time  t=0.23 


Fig.  9.  Plot  of  average  vorticity  magnitude  of  the  thirty  control  points  with  highest 
vorticity  magnitude  for  cases  with  strong,  moderate  and  weak  vortices,  with 
r  /  Ina^U  =  1.6  (dashed  curve),  0.79  (dashed-dotted  curve)  and  0.32  (solid 

curve),  in  comparison  to  RDT  results  (circles).  The  RDT  results  are  the  same 
for  all  values  of  F  /  2noM. 


The  average  vorticity  magnitude,  of  the  30  control  points  with  the  highest 
vorticity  magnitudes,  divided  by  the  initial  maximum  vorticity  magnitude  is  plotted 
in  figure  9  versus  time  for  all  three  values  of  the  parameter  T  UtiCqU  with  the  high 
resolution  runs.  The  time  interval  shown  corresponds  to  that  over  which  the  vortex 
penetration  occurs.  The  selection  of  30  control  points  over  which  to  average  vorticity 
magnitude  was  found  to  give  values  fairly  close  to  the  peak  vorticity  values,  while  being 
somewhat  smoother  than  averages  with  smaller  number  of  points.  Results  for  the  full 
inviscid  calculations  are  shown  by  solid,  dashed-dotted,  and  dashed  curves  for  cases  with 
r  /  ItiUCq  =0.32,  0.79  and  1.6,  respectively.  The  results  for  /  co,„„  obtained  from 

the  three  RDT  calculations  were  identical,  since  the  vorticity  was  deformed  only  by  the 
ambient  potential  flow.  The  RDT  results  are  shown  in  figure  9  by  circular  symbols.  The 
values  of  /  (j0,„„  shown  in  figure  9  for  the  two  full  inviscid  calculations  with 

circulation  of  F  =  0.2  and  0.5  are  very  close  for  the  RDT  results.  In  the  case  with  F  =  1 , 
the  value  of  this  ratio  deviates  from  that  for  the  RDT  calculation  at  a  time  of  about  0.31, 
after  about  80%  of  the  blade  penetration  is  complete.  The  deviation  of  to^  /  0),„„  for  the 

strong  vortex  case  from  that  obtained  for  weaker  vortices  could  be  caused  by  stretching  of 
the  vorticity  by  interaction  with  its  image  in  the  blade,  which  would  not  be  as  prevalent 
for  cases  with  weaker  vortices,  or  by  the  numerical  difficulties  encountered  at  the  end  of 
this  run,  as  noted  in  the  discussion  of  figure  7.  It  is  also  noted  that  the  leveling  off  of  the 
maximum  vorticity  values  in  figure  9  after  the  penetration  is  about  75%  complete  (for 
t  >  0.30)  is  due  to  the  limitation  on  the  ability  of  the  vorticity  blobs  to  resolve  the  very 
thin  vortex  ribbon  about  the  blade  surface  during  the  latter  part  of  the  runs. 

(in)  Blade  surface  pressure  and  force 

A  contour  plot  of  the  change  in  blade  surface  pressure  from  its  value  without  the 
vortex  present  is  given  in  figure  10  for  a  moderate  resolution  run  with  F  /  IncfU  =  1.6. 

The  plot  is  made  at  t  =  0.26,  at  which  time  the  vortex  penetration  is  roughly  half 
completed,  and  shows  the  vorticity  vector  only  at  every  fourth  control  point  in  order  to 
indicate  the  general  location  of  the  vortex.  The  pressure  signature  on  the  blade  surface  is 
dominated  by  a  'V  of  low  pressure,  in  the  center  of  which  (at  the  blade  leading  edge) 
there  is  a  patch  of  high  pressure.  The  low  pressure  patches  on  the  top  and  bottom  of  the 
blade  are  due  to  the  low  pressure  within  the  vortex  core,  and  the  low  and  high  pressure 
regions  along  the  blade  leading  edge  are  caused  by  the  induced  velocity  from  the  vortex 
as  it  impinges  on  the  blade.  The  "waviness"  of  the  element  arrows  at  the  top  and  bottom 
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Fig.  10.  Contour  plot  of  the  pressure  change  on  the  blade  surface  due  to  the  vortex  for 
a  case  with  F  /  2kOqU  =  1.6  and  at  time  t  =  0.26.  One-fourth  of  the  vorticity 

control  point  positions  are  shown  to  indicate  the  position  of  the  vortex. 
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of  the  ring  is  an  optical  effect  due  to  looking  along  the  core  axis  (such  waviness  is  not 
apparent  in  the  side  view  in  figure  8). 


The  pressure  signature  for  the  weak  vortex  case  appears  similar  to  that  shown  in 
figure  10,  except  that  high  pressure  regions  are  also  found  on  the  top  and  bottom  of  the 
blade  on  the  opposite  side  (i.e.,  for  opposite  values  of  z)  of  that  on  which  the  low  pressure 
regions  occur.  The  general  appearance  of  the  pressure  signature  for  weak  vortices  is  that 
of  two  interlocking  Vs  of  high  and  low  pressure. 

A  plot  of  time  variation  of  the  pressure  coefficient  associated  with  the 

minimum  surface  pressure  on  the  blade,  defined  by 


C„„.„  = 


p,Tnin 


Pq  Pmin 


(1) 


is  given  in  figure  1 1  for  the  high  resolution  calculations  with  three  different  values  of  ring 
circulation.  Results  of  the  full  inviscid  calculations  are  shown  by  continuous  curves  and 
results  for  RDT  calculations  are  shown  by  symbols.  In  the  implementation  of  RDT  in  our 
computations,  the  vorticity  is  convected  with  only  the  ambient  potential  flow  but  the  full 
velocity  is  used  for  evaluation  of  pressure.  It  is  found  that  if  only  the  ambient  potential 
flow  is  used  in  the  pressure  calculation,  the  computed  unsteady  drag  force  on  the  blade  is 
much  too  small.  In  an  analytical  study  using  RDT  together  with  the  thin  airfoil 
approximation,  Howe  (1989)  also  found  that  the  spanwise  velocity  components,  induced 
by  the  vortex  and  its  image  over  the  blade  surface,  play  an  important  role  in  determination 
of  the  unsteady  blade  force. 

The  results  of  the  full  inviscid  computations  for  the  lower  two  circulation  values 
compare  closely  with  the  RDT  runs,  but  the  results  for  the  highest  circulation  value 
oscillate  about  the  RDT  results.  The  pressure  field  for  runs  with  strong  vortices  is 
dominated  by  the  vortex-induced  velocity  and  would  therefore  be  expected  to  vary 
approximately  in  proportion  to  T^.  For  much  weaker  vortices,  the  pressure  field  is 
influenced  more  by  the  perturbation  which  the  vortex  causes  to  the  ambient  velocity  field, 
and  would  therefore  vary  approximately  in  proportion  to  UT.  Since  the  three  cases 
plotted  in  figure  1 1  include  vortices  ranging  from  fairly  weak  to  fairly  strong  (compared 
to  the  free  stream  flow),  the  minimum  pressure  coefficient  is  lower  for  cases  with 
stronger  vortices. 
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Fig.  11.  Plot  of  variation  of  the  minimum  value  of  the  blade  surface  pressure  coefficient 
Cp  min  with  time  for  cases  with  strong,  moderate  and  weak  vortices, 

r  /  2%GqU  =  1.6  (dashed  curve),  0.79  (dashed-dotted  curve)  and  0.32  (solid 

curve),  in  comparison  to  RDT  results,  denoted  by  circles,  triangles  and  squares, 
respectively. 
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The  computational  results  for  vortex-induced  drag  D  on  the  blade  are  plotted  in 
figure  12  for  high-resolution  runs  in  terms  of  the  drag  coefficient  C^,  defined  by 


-D 


(2) 


The  computed  results  for  from  the  RDT  calculations  are  found  to  be  identical  (to 
within  plotting  accuracy)  for  all  three  values  of  circulation  considered,  even  though  the 
pressure  results  shown  in  figure  1 1  are  quite  different.  Furthermore,  the  results  of  the  full 
in  viscid  calculations  are  found  to  coincide  very  closely  to  the  RDT  results  for  t  <  0.29 
for  all  three  circulation  values  considered.  For  f  >  0.29,  the  full  in  viscid  calculations  for 
the  cases  with  F  =  0.5  and  1.0  (corresponding  to  e  =  0.79  and  1.6,  respectively)  give 
somewhat  larger  values  for  Cp  than  predicted  by  RDT.  Since  both  the  maximum 
vorticity  values  and  the  minimum  pressure  values  for  the  case  with  T  =  0.5  agree  closely 
with  the  RDT  results  during  this  time  period  (as  shown  in  figures  9  and  11),  it  seems 
unlikely  that  the  observed  differences  in  drag  coefficient  for  this  case  are  due  to 
numerical  resolution  difficulties  (which  would  tend  to  have  most  effect  on  points  with 
largest  vorticity  values  in  the  ribbon  about  the  blade  surface).  These  differences  in  Cp 
with  the  RDT  predictions  seem  to  arise  from  breakdown  of  the  RDT  assumption  as  the 
parameter  e  approaches  and  exceeds  unity,  resulting  in  differences  in  large-scale 
deformation  of  the  vortex  between  the  RDT  and  full  inviscid  calculations. 


The  variation  of  blade  penetration  force  with  T^,  which  is  evident  from  figure  12 
during  the  first  70-75%  of  the  penetration  process,  can  be  explained  by  recalling  the 
identity  obtained  by  Howe  (1989,  Eq.  (A.4)),  which  gives  the  drag  force  due  to  an 
inviscid  vorticity  field  external  to  a  stationary  body,  with  constant  uniform  flow  at 
infinity,  by  the  integral 


D  =  pJ  V(^-(roxu,„^)dv.  (3) 

V 

Here  ([)  is  the  velocity  potential  due  to  a  uniform  potential  flow  (with  unity  flow  speed  at 
infinity)  past  the  body  and  is  the  velocity  induced  by  the  vortex  and  its  image.  If  the 

deformation  of  the  vorticity  field  were  due  primarily  to  the  ambient  potential  flow  past 
the  body,  as  assumed  in  RDT,  then  both  c5  and  would  be  proportional  to  F,  so  that  D 
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Fig.  12.  Plot  of  variation  of  the  drag  coefficient  with  time  for  cases  with  strong, 

moderate  and  weak  vortices,  F  /  2kCqU  =  1.6  (dashed  curve),  0.79  (dashed- 

dotted  curve)  and  0.32  (solid  curve),  in  comparison  with  RDT  results  (circles). 
The  RDT  results  are  the  same  for  all  values  of  F  /  ItiCqU. 
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would  be  proportional  to  and  (3)  would  yield  the  result  that  the  value  of  is 
independent  of  F  (as  our  computed  results  with  RDT  show).  On  the  other  hand,  if  the 
vortex  were  sufficiently  strong  that  the  vortex-induced  or  image  velocity  contributed 
significantly  to  deformation  of  the  vorticity  field,  then  the  vortex  stretching  term  in  the 
vorticity  transport  equation  would  exhibit  both  terms  which  are  linear  and  quadratic  in  F, 
and  (3)  would  indicate  that  the  drag  will  no  longer  be  proportional  to  F^. 

(iv)  Comparison  of  different  vortex-induced  forces  on  the  blade 

It  may  be  of  interest  at  this  point  to  compare  the  magnitude  of  the  blade 
penetration  force  with  that  of  other  forces  present  in  typical  normal  BVI  applications.  For 
instance,  for  non-zero  angles  of  attack  a,  the  normal  BVI  force  is  calculated  in  Marshall 
and  Yalamanchili  (1994)  prior  to  impact  of  the  blade  on  the  vortex  using  a  vortex 
filament  model.  This  force  is  found  to  increase  with  increasing  values  of  a  and  T  !  <3 q. 

For  instance,  in  figure  10  of  Marshall  and  Yalamanchili  (1994)  (which  is  for  a  case  with 
T/cTo  =  5  and  attack  angles  of  a  =  0°  and  15°),  the  peak  drag  coefficients  are  found  to  be 

about  =  0.3  and  0.5,  respectively.  In  comparison,  if  we  assume  that  the  blade 
penetration  force  increases  in  proportion  to  the  ratio  T !  (5 q  (for  sufficiently  small  values 
of  this  ratio)  and  vanishes  for  F  /  cJq  =  0  (both  of  which  are  consistent  with  the  result  in 

Eq.  (3.20)  of  Howe,  1989),  the  peak  value  from  figure  12  can  be  extrapolated  for  a  case 
with  T/rto  =  5  to  give  a  peak  value  of  approximately  Co=0.5.  For  cases  with 
F/ao  =  0(1),  the  blade  penetration  force  thus  seems  to  make  a  significant  contribution 
to  the  total  vortex-induced  force  on  the  blade. 


3.  Conclusions  and  Future  Work 


Three  regimes  of  vortex-blade  interaction  have  been  experimentally  identified, 
based  on  the  values  of  the  two  dimensionless  parameters  F  /  a  and  F  /  IncsU .  For  large 
values  of  F  HtvsU  ,  the  interaction  is  controlled  by  vorticity  eruption  from  the  blade 
boundary  layer,  followed  by  wrapping  of  the  secondary  vorticity  about  the  primary  vortex 
and  subsequent  breakup  of  the  primary  vortex.  Calculations  of  vortex-induced  separation 
from  a  circular  cylinder  have  been  recently  been  reported  by  Affes  et  al.  (1994)  based  on 
boundary  layer  theory.  However,  there  is  very  little  literature  pertaining  to  the  problem  of 
wrapping  of  secondary  vorticity  about  the  primary  vortex  and  the  observed  primary 
vortex  breakup  in  this  regime. 
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As  T  UtigU  decreases  below  about  10-20,  the  secondary  vorticity  plays  a 
progressively  smaller  role  in  the  vortex-blade  interaction.  For  the  case  of  small 
T UtioU  and  7/ a  much  greater  than  0(1),  the  vortex  is  observed  to  bend  about  the 
blade  (or  cylinder)  much  as  is  predicted  by  the  vortex  filament  theory.  When  vortex 
breakup  finally  does  occur  in  this  regime,  it  is  preceded  by  formation  of  short  helical 
waves  on  the  vortex  core,  which  grow  and  impinge  on  the  body  surface.  If  these  waves 
are  due  to  inviscid  instability  of  the  primary  vortex  with  its  image,  they  should  be 
observed  in  the  vortex  filament  predictions  provided  that  the  spatial  resolution  is 
sufficient  and  that  initial  small  perturbations  over  a  spectrum  of  wavelengths  are 
introduced  on  the  vortex  core.  Alternatively,  the  waves  could  be  due  to  interaction  with 
shed  vorticity  from  the  blade  boundary  layer  as  the  primary  vortex  is  bent  about  the  body. 
Based  on  estimates  from  linear  stability  theory  and  flow  visualization,  it  seems  that  both 
of  these  effects  could  be  significant. 

For  the  case  of  small.  r/27CCTC/  and  7/ a  of  0(1)  or  smaller,  the  blade 
penetrates  into  the  vortex  core.  The  deformation  of  the  vorticity  field  during  blade 
penetration  was  examined  in  detail,  for  an  inviscid  flow,  by  Marshall  and  Grant  (1995). 
The  vorticity  is  observed  to  accumulate  into  a  thin  vorticity  "ribbon"  which  wraps  about 
the  blade  surface.  Comparison  of  results  of  direct  computations  from  the  Euler  equations 
and  calculations  using  the  RDT  approximation  indicates  that  RDT  is  valid  for  this 
problem  for  (F  /  2koU){T  /a)  less  than  about  unity.  For  values  7  /  o  much  less  than 
unity,  the  blade  penetration  force  would  be  small  compared  to  other  forces  present. 
Problems  which  have  not  been  resolved  for  this  regime  include  the  decay  of  the  vorticity 
"ribbon"  due  to  viscous  effects  and  the  force  exerted  by  blockage  of  the  axial  flow  within 
the  vortex  as  it  is  cut  by  the  blade. 
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